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THE  MANN  GRUBBS-METHOD 

Nancy  R.  Mann* 
Frank  E.  Grubbs 


or _ 

n1 5-.tr  i  tut  1  ,a ! 


«v  ti  fc».» 


Consider  a  series  system  made  up  of  k  independent  subsystems, 
each  having  exponentially  distributed  failure  time  T  with 
either  no  censoring  of  the  subsystem  data  or  only  type-II 
censoring.  Testing  is  without  replacement.  For  this  model, 
system  reliability  Rs (t^)  at  time  t^  >  0  is  equal  to 


II  exp (-t  A.)  -  exp(-t  EX), 


i=l 


m  j 


°j=l  j‘ 


(1) 


with  Aj  >  0,  j=l,...,k,  the  hazard  rate  for  the  jth  sub¬ 
system.  One  can  determine  an  upper  confidence  bound  on  the 

k 

system  hazard  rate  $  =  I  A.,  and  from  this,  of  course,  can 

jWl  J 

then  obtain  a  lower  confidence  bound  on  R  (t  )  «  exp(-<}>t  ). 

s  m  m 

The  Mann-Grubbs  (1972)  method  for  the  determination  of 
confidence  bounds  for  system  reliability  has  been  adapted  to 
several  situations,  but  was  derived  originally  for  this  par¬ 
ticular  model.  For  this  exponential  series-system  model  with 
type-II  censoring,  the  method  yields  a  lower  confidence  bound 
that  very  closely  approximates  the  lower  confidence  bound  that 
is  most  accurate  (has  the  highest  probability  of  being  close 
to  the  true  system  reliability)  for  all  values  of  system  re¬ 
liability,  among  exact  confidence  bounds  that  are  unbiased. 

The  restriction  of  unbiasedness  is  necessary  here  because  of 
the  nuisance  parameters  A^,***!^,  the  hazard  rates  for  the  k 
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independent  subsystems.  For  binomial  models  and  exponential 
models  with  censoring  by  time,  the  Mann-Grubbs  (1974)  method 
uses  an  approach  similar  to  that  used  for  the  exponential 

j 

model  with  type-11  censoring.  However,  because  of  discrete¬ 
ness  of  the  data  in  the  binomial  case,  optimal  bounds  are  1 
more  difficult  to  approximate.  For  all  of  these  situations, 

J 

however,  the  method  is  quite  simple  to  implement. 

For  the  model  given  by  equation  (1)  with  k  =  2,  Lentner 
and  Buehler  (1963)  derived  the  uniformly  most  accurate  un¬ 
biased  lower  confidence  bound  on  R  (t  ) .  Generalization  to 

s  m 

k  >_  2  was  made  by  El  Mawaziny  (1965)  in  his  doctoral  thesis. 

As  noted  above,  their  results  depend  upon  the  assumption  that 

for  the  jth  subsystem  ^prototypes  have  been  tested  until 

r^,  with  1  Tj  Jl  n j  >  failures  occur,  j  =  1 . k.  For  the 

jth  subsystem,  one  observes  the  ith  smallest  failure  times 

t.  ,  i=l,...,r..  One  then  computes  the  total  time  on  test 
J  J  £ 

for  the  jth  subsystem,  w  «■  2^  t  +  (n.-r  )t  ,  j=l . 

J  i=i  J  J  rj»J 

k.  Calculation  of  a  lower  confidence  on  R  (t  )  based  on  the 

S  IQ 

w, 's  and  the  r  *s  by  El  Mawaziny's  method  must  be  performed 

J  J 

iteratively  by  means  of  a  computer.  If  both  the  number  of 
subsystems  and  the  total  number  of  failures  are  large,  prob¬ 
lems  of  loss  of  precision  will  result.  See  Mann  (1970). 


The  Mann-Grubbs  (M-G)  method  eliminates  the  need  for  a 
computer  and  the  problem  of  loss  of  precision  resulting  from 
the  lengthy  calculations.  The  approach  is  based  on  the  fact. 
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demonstrated  by  Mann  and  Grubbs  (1972),  that  the  posterior 

k 

distribution  of  $  =  £  X  is  that  of  a  sum  of  weighted  non- 

j=l  J 

central  chi  squares.  This  sum  can  be  well  approximated  for 

present  purposes  by  a  single  weighted  chi-square  variate  with 

mean  m  and  variance  v.  Thus,  one  assumes,  using  a  two-moment 

fit  to  chi-square,  that  2m$/v  is  a  chi-square  variate  with 

2m2/v  degrees  of  freedom  (see  Patnaik,  1949). 

The  expressions  for  the  conditional  mean  m  and  variance  v 

of  the  system  hazard  rate  derived  by  Mann  and  Grubbs  for  this 

model  have  been  simplified  by  Mann  (1974)  and  are  given  by 
k  (r.-l)  , 


j-1  Wj 


k  (t.-l) 

I  - 

j-l  Wj 


where  w^  is  the  smallest  of  the  w^'s. 

Once  the  hazard-function  moments  have  been  calculated, 
then  the  Wilson-Hilferty  transformation  of  chi-square  to 
normality  can  be  used  to  facilitate  the  calculations,  since 
the  number  of  degrees  of  freedom  v  =  2m2/v  for  the  approxi¬ 
mate  chi-square  variate  2m $/v  is  not  generally  an  integer. 
To  approximate  the  uniformly  most  accurate  lower  confidence 


bound  R  (t  )  on  series-system  reliability  R(t  )  at  time  t  and 
s  m  mm 

at  confidence  level  1  -  a  (incorporating  the  Wilson-Hilferty 
transformation),  one  calculates 
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R  (t  )  =  exp[-t  m(l  -  — + 


is 

z.  v  3 


where  z^  is  the  lOOyth  percentile  of  a  standard  normal  dis¬ 
tribution.  The  Wilson-Hilf erty  transformation  yields  an 
approximation  to  chi  square  in  this  context  that,  for  three 
or  more  degrees  of  freedom,  is  accurate  to  within  a  unit  or 
two  in  the  second  significant  figure. 

For  an  example  of  calculation  of  an  approximate  confi¬ 
dence  bound  R  (t  )  on  series  system  reliability,  we  consider 
— s  m 

an  independent  series  system  containing  three  subsystems, 

each  with  exponential  failure  time.  For  each  of  the  three 

subsystems,  prototypes  have  been  life  tested  with  type-II 

censoring  resulting  in  total-times-on  test  w1  =  42.753, 

w2  =  45.791  and  w3  =  31.890,  with  r3  =  4,  r2  •=  3,  and 

r3  =  2.  To  obtain  R^t^)  using  (2),  (3)  and  (4),  one  forms 

4-1  3-1  .  2-1  1 

m  42.753  45.791  31.890  31.890  " 


v  =  ftzl  +  J-i-  +  ±Z±  +  i 

(42.753)2  (45.7  91) 2  (31.890)2  (31.890)2 

=  .00456. 


Thus,  an  approximate  ninety  percent  lower  confidence  bound 

at  t  B  1  is  calculated  as 
m 

R  ( 1) -exp { - .  17 656[  1- -  — -{’-17-5- 2  +  1.282/.00456/.52968]3} 

— s  “ 

■  .766  . 

The  F.l  Mawaziny  optimal  exact  lower  confidence  bound  on 
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R(l)  is  .772.  See  Mann  and  Grubbs  (1972)  or  Mann,  Schafer, 
s 

Singpurwalla  (1974). 

Schoenstadt  (1980)  uses  simulation  to  compare  an  exact 
procedure  of  Lieberman  and  Ross  (1971)  and  this  approximate 
method  for  obtaining  lower  confidence  bounds  on  system  reli¬ 
ability  under  the  model  specified  by  equation  (1).  He  con¬ 
cludes  that  the  results  of  the  simulation  runs  "seem  to 
demonstrate  the  superiority  of  the  M-G  bounds  in  those  in¬ 
stances  that  can  be  considered  of  practical  importance,"  and 
also  demonstrates  that  the  M-G  bounds  are  exact  to  the 
accuracy  of  his  simulation. 

For  binomial  data,  the  parameter  of  interest  for  which  a 
posterior  mean  and  variance  (conditional  on  the  failure  data) 
are  calculated  is  £,  the  negative  of  the  logarithm  of  system 
reliability  R  or,  for  parallel  systems,  the  probability  of 
system  failure,  1  -  R.  The  parameter  £  has  been  demonstrated 
by  Mann  (1974a)  to  have  a  posterior  distribution  that  is 
approximately  proportional  to  a  chi-square  variate,  as  does 
$  in  the  exponential  model.  The  expressions  for  the  condi¬ 
tional  mean  and  variance  of  £  resemble  the  expressions  (2) 
and  (3)  only  for  models  for  which  randomized,  bounds,  commonly 
used  in  binomial  models,  are  obtained. 

Expressions  for  obtaining  approximate  lower  confidence 

bounds  on  R  or  R  (either  randomized  or  nonrandom  zed)  for 
s  p 

the  binomial  model  and  for  an  exponential  model  with  fixed 
censoring  times  can  be  found  in  Section  10.4  Mann,  Schafer, 
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Singpurwalla  (1974).  Comparisons  of  results  with  -hose  of 
other  methods  are  also  given  this  for  all  applicable  models, 
as  well  as  methods  for  combining  results  to  obtain  bounds  fo'r 
models  more  general  than  simple  series  or  parallel  systems. 
The  latter  are  also  discussed  in  Mann  and  Grubbs  (1974). 
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